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We study the wetting transition and the directed polymer derealization transition on diamond hi- 
erarchical lattices. These two phase transitions with frozen disorder correspond to the critical points 
■ of quadratic renormalizations of the partition function. ( These exact renormalizations on diamond 

lattices can also be considered as approximate Migdal-Kadanoff renormalizations for hypercubic 
f^**; , lattices). In terms of the rescaled partition function z = Z/Ztyp, we find that the critical point 

. corresponds to a fixed point distribution with a power-law tail P c (z) ~ $(ln z)/z 1+)1 as z — > +00 ( 

up to some sub-leading logarithmic correction $(lnz)), so that all moments z™ with n > fj, diverge. 
, For the wetting transition, the first moment diverges ~z = +00 (case < fj, < 1), and the critical 

temperature is strictly below the annealed temperature T c < T ann . For the directed polymer case, 
the second moment diverges z 2 — +00 (case 1 < fj, < 2), and the critical temperature is strictly 
below the exactly known transition temperature T2 of the second moment. We then consider the 
correlation length exponent v : the linearized renormalization around the fixed point distribution 
coincides with the transfer matrix describing a directed polymer on the Cayley tree, but the random 
weights determined by the fixed point distribution P c {z) are broadly distributed. This induces some 
I , changes in the travelling wave solutions with respect to the usual case of more narrow distributions. 
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I. INTRODUCTION 

A. Real-space renormalizations for disordered systems 
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The choice to work in real-space to define renormalization procedures, which already present a great interest for 
pure systems [l|, becomes the unique choice for disordered systems if one wishes to describe spatial heterogeneities. 
Whenever these disorder heterogeneities play a dominant role over thermal or quantum fluctuations, the most ap- 
propriate renormalizations are strong disorder renormalizations Q introduced by Ma-Dasgupta Q : as shown by 
t-H ■ Fisher 0], these strong disorder renormalization rules lead to asymptotic exact results if the broadness of the disor- 
der distribution grows indefinitely at large scales. However, for disordered systems governed by finite-disorder fixed 
points, where disorder fluctuations remain of the same order of thermal fluctuations, one needs to use more standard 
real-space renormalization procedures, such as Migdal-Kadanoff block renormalizations They can be considered 
in two ways, either as approximate renormalization procedures on hypercubic lattices, or as exact renormalization 
• ■ procedures on certain hierarchical lattices @, 0] • One of the most studied hierarchical lattice is the diamond lattice 
which is constructed recursively from a single link called here generation n = (see Figure [I} : generation n = 1 
consists of b branches, each branch containing 2 bonds in series ; generation n = 2 is obtained by applying the same 
transformation to each bond of the generation n = 1. At generation n, the length L n between the two extreme sites 



A and B is L„ = 2™, the total number B n of bonds is B n = (26)™ = Ln" ff so that d e ff(b) = 1 ", 2 ^ represents some 



In 2 

. !Zh ! effective dimensionality. 

On this diamond lattice, various disordered spin models have been studied, such as for instance the diluted Ising 
model Q , random bond Potts model , and spin-glasse s [Toll . Disordered polymer models have also been considered, 
in par ticular the wetting on a disordered substrate [HI, Il2|| and the directed polymer model [H, [13, [H, [H, [U EH 
[l9l. l2fj. l2l| . In this paper, we focus on these two polymer models that are described by quadratic renormalization of 
their partition functions as we now recall. 

B. Wetting transition with disorder on the diamond lattice 

On the diamond lattice, the adsorption of a polymer on a disordered substrate is described by the following quadratic 
recursion for the partition function Z n of generation n [Tlj 

Z n+1 = Z^zW + (b-l)Y* (1) 
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FIG. 1: Hierarchical construction of the diamond lattice of branching ratio 6. 



where Y n = b Ln 1 represents the number of walks between the two extreme points and satisfies the recursion without 
disorder 

Y n +i = bY* (2) 

and where Zn^ and Zn' represent two independent copies of generation n. At generation n — 0, the lattice reduces 
to a single bond with a random energy e, for instance drawn from the Gaussian distribution 

*) - 4-"* (3) 



and thus the initial condition for the recursion of Eq. [T] is simply 

Z n=0 = e-P* (4) 
The temperature only appears in this initial condition. 



C. Directed polymer on the diamond lattice 



The model of a directed polymer in a random medium [22j can also be studied on the diamond hierarchical lattice 
EH El, El, E3, El, El, El, EH . The partition function Z n of the n-generation satisfies the exact recursion [3] 

Z n+l =j^Z^Z^ (5) 

a=l 

where (Z^ , zi 2b) ) are (26) independent partition functions of generation n. At generation n — 0, the lattice 
reduces to a single bond with a random energy e, for instance drawn from the Gaussian of Eq. [3] and thus the initial 
condition for the recursion of Eq. [5] is again given by Eq. 01 



D. Organization of the paper 



In this paper, we study the critical points of the quadratic renormalizations described above that correspond 
to derealization transitions for the polymer. These two transitions are of course different in nature, since the 
wetting transition already exists in the pure case, whereas the directed polymer transition only exists in the presence 
of disorder. However, we will show below that the quadratic form of the renormalizations induce some common 
properties. It is thus interesting to study them along the same lines to stress their similarities and differences. The 
paper is organized as follows. The wetting transition on a disordered substrate is discussed in Section m and studied 
numerically in Section IIIIl The directed polymer transition is discussed in Section IIV1 and studied numerically in 
Section [V] In Section IVTl we compare the results on the diamond lattice with respect to the same disordered polymer 
models defined on hypercubic lattices. Section [VIII contains the conclusion. The Appendix [Al contains a reminder on 
multiplicative stochastic processes which is used in Sections [TT] and IIVI 
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II. WETTING ON A DISORDERED SUBSTRATE 



To study the wetting recursion, it is convenient to introduce the reduced partition function z n and the associated 
free- energy /„ defined by [HI 

zn = ^ - e"". (6) 



to rewrite the recursion of Eq. [T] as 

z n+l = 7 (') 



A. Reminder on the pure case 

In the pure case, the ratios z n defined in Eq. [5] are not random but take a single value R n , and the recursion of Eq. 
[7] reduces to a one-dimensional mapping T 

R n+ i = Rl+ ^ 1 = T(R n ) (8) 

discussed in [ll| : for b > 2, there exists two attractive fixed points Roc = 1 (delocalized phase) and Roo = +oo 
(localized phase) separated by the repulsive fixed point R c (critical point) with 

R c = b-l (9) 

The critical exponents are determined by the linearization of the recurrence around the fixed point. Setting 
R n = R c + S n , one obtains at linear order 

<y„+i - with A = T'(R C ) = ^ = ~ (10) 

Note that this factor A = T'{R C ) > 1 describing the instability of the critical point also governs the growing of the 
energy E n exactly at criticality [201 ] . since the recursion for the energy 

Rl(2E n ) 



becomes at criticality 



E n+ i(T c ) = R2 i E n (T c ) = XE n (T c ) (12) 



To understand why the same factor A appears, one may introduce the product U n = R n E n that satisfies the recursion 

Un+1 = 2 -^U n (13) 
o 

It is then clear that at criticality it coincides with the recursion of the variables S n (Eq. [T0|) 

In conclusion, the varia 
with the critical exponent 



In conclusion, the variable 8 n or the energy E n at criticality grows as A™ = L\[ v in terms of the length L n = 2 r ' 



_ In 2 _ In 2 

V ~ hiA ~ ln T'{R C ) (M) 

The specific heat exponent satisfies the hyperscaling relation a = 2 — v. We refer to [ll| for more details. 
Let us now summarize the changes that the presence of frozen disorder will induce : 

(i) the one-dimensional mapping of the pure case R n +i — T(R n ) will become the iteration of a probability distri- 
bution Q n+ x(z) = T{Q n (z)} 

(ii) the critical value R c of the pure case will become an invariant probability distribution Q c (z) — T{Q c {z)} 

(iii) the critical exponent v determined by the derivative T'(R C ) in the pure case will be determined by the linearized 
iteration around the fixed point distribution Q c (z). 

But before concentrating on the critical point, we first describe the properties of the renormalization group (RG) 
flow with disorder in the limits of high and low temperatures. 
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B. High-temperature RG flow 

In the high temperature phase, the variables z n defined in Eq. [6] flow towards 1 or equivalently the free-energies /„ 
decay to zero. The linearization of the recursion in this regime yields 

f (1) , f (2) 

Jn+1- £ (15) 

For b > 2, this high-temperature phase exists, the probability distribution of the free-energy converges to a Gaussian, 
the average and the width decays as power-laws of the length L n = 2™ 

_ in 4 

fn «L„'» ! (16) 



/•?-(/„) 2 «L r T^ (b) with co' w (b) = ^ (17) 

C. Low-temperature RG flow 

In the low-temperature phase, the free-energies /„ of Eq. [5] grow extensively with the length L n — 2", and thus at 
large scale, the recursion is dominated by the first term in Eq. [7] 

/„ +1 ^/W+/( 2 > + ... (18) 

The probability distribution of the free-energy thus converges to a Gaussian, the average and the width grows as 

Jn^L n (19) 



/*-(/n) 2 0Ci# a (20) 

D. Analysis of the critical invariant distribution 

At criticality, to avoid the high-temperature and low-temperature described above, the free-energy /„ of Eq. [6] 
should remain a random variable of order 0(1) with some scale-invariant probability distribution P c (f) defined on 
] — oo,0]. Equivalently, the variable z n = e~^ c ^ n should have a scale-invariant probability distribution Q c (z) defined 
on [1, +oo [. In the following, we derive some of their properties. 

1. Left-tail behavior of the free-energy distribution 

Let us introduce the left-tail exponent r\ c 

lnP c (/) ~ - 7 (-/)n. + ... (21) 

f—>—oo 

where (...) denote the subleading terms. 

In the region where / — > — oo, one has effectively the low-temperature recursion 

/ ~ / « +/ (2) + ... (22) 

A saddle-point analysis shows that if and have a probability distribution with the left tail given by Eql2T| 
their sum / has for left tail lnP c (/) ~ — j(— /) ? ' c 2 1_, ' c + .... The stability of the critical distribution thus fixes the 
value of the left tail exponent to 

Vc = 1 (23) 

So the distribution P c {f) decays exponentially 

P c (f) ~ etf(...) (24) 

/— >— oo 
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This means that the corresponding distribution Q c (z) of z = e presents a power-law tail 

$(hiz) 

with some exponent 



M = ^ (26) 



and where $(lnz) represents the subleading terms. 



2. Analysis in terms of multiplicative stochastic processes 

The fact that a power-law appears in the stationary distribution of some random iteration is reminiscent of mul- 
tiplicative stochastic processes, whose main properties are recalled in Appendix [A] For a multiplicative stochastic 
process X n described by Eq. IA11 the stationary distribution presents a power-law tail of exponent [i that can be 
computed in terms of the statistics of the random coefficient a n via Eq. IA41 Here, for the quadratic renormalization of 
Eq. it is the process z n itself that also plays the role of the random multiplicative coefficient. As a consequence, it 
is instructive to analyse the recursion of Eq. [Jj along the same lines used to study multiplicative stochastic processes. 

The necessary stability condition of Eq. IA2I translates here into the following condition 



ln-=/ dzQJz)lnz-lnb< (27) 
b Ji 

The condition of Eq. IA4I that ensures the stability of the power-law tail via iteration translates into the following 
self-consistent condition for the tail exponent jj. introduced in Eq. [23 



2(|)" =2J dzQ c (z) (£f = l (28) 

The argument is similar to the computation of Eqs IA5HA61 the additional factor of 2 coming from the fact that z 
large corresponds to either z^ 1 ' large or z^ large. The condition of Eq. [28] means in particular that the subleading 
term $(lnz) in Eq. 1251 should ensure the convergence at (+oo) of the following integral 

z» = / dzQ c {z)z tl ~ / — $(ln2r)~ / dw<S>{w) < +oo (29) 
Since $ is a subleading term in Eq 1251 it should not contain an exponential, so its decay should be a power-law 

$(w) ~ — r— with a > (30) 

v ' w-> oc W 1+a 

Then moments of order k < fi are finite, whereas moments of order k > \l diverge 

dzQ c (z)z k = +oo for k > fi (31) 



In contrast with multiplicative stochastic processes where the condition of Eq. IA4I allows to compute the tail 
exponent ji in terms of the known statistics of the random coefficient a n , we have obtained here only a self-consistent 
equation : the selected tail exponent /j, in the region z — > oo is the exponent that satisfies the condition of Eq. [^5] 
that involves the whole distribution for z € (1, +oo[. However, even if we cannot explicitly compute this exponent /x, 
we can try to locate it with respect to integer values by considering the integer moments. 

3. Reminder on transitions of integer moments 

An important property of quadratic renormalizations is that they lead to closed renormalizations for the integer 
moments. We now briefly recall the behavior of the first moments discussed in (Tlj . The closed recursion satisfied by 
the first moment [ll| 

___ _ (~z^) 2 + ^ - 1 (32) 
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coincides with the pure case equation of Eq. [8J Using the initial condition of Eq. [4] the unstable fixed point of Eq. 
fallows the define the annealed temperature via e~ ei / T °- n ™ =6—1 : for T > T ann , the averaged value goes to 1, 
whereas for T < T ann , the averaged value goes to +oo. So T ann represents the transition of the first moment. To 
locate T c with respect to T ann , we have to distinguish two possibilities 

(i) if the tail exponent fj. of Eq. [^satisfies < fi < 1, then its first moment diverges at criticality = +oo and 
we have the strict inequality T c < T ann . 

(ii) if the tail exponent /x satisfies fj, > 1, then its first moment is finite at criticality. The only possible finite stable 
value is z^ = b — 1 and the critical temperature then coincides with the annealed temperature T c = T ann . However, 
the analysis of the recursion for the variance leads to the conclusion that the critical temperature is strictly lower 
than the annealed temperature T c < T ann as soon as disorder is relevant b > 2 + y/2 ~ 3.414 [Tlj . 

In conclusion, whenever disorder is relevant at criticality, one has the strict inequality T c < T ann , the first moment 
diverges z c = +oo, and the tail exponent fi of Eq. [25] is smaller than 1 

< n < 1 (33) 



E. Critical exponent v 

1. Equivalence with a directed polymer on a Cayley tree 

In the pure case, the critical exponents are obtained from the linearization around the fixed point (see section III Ajl . 
To follow the same strategy in the disordered case, we set z n — z c + S n . At linear order, we obtain the recursion 

(1) (2) 

S n + 1 ~ ^) + zc_ 5{ i) (34) 

(1 2) 

where z c ' are distributed with the critical distribution Q c (z). As in the pure case, it is also interesting to write the 
recursion for the energy E n 

En+1 ~ £>£> + h-i (35) 

so that the combination U n = z n E n satisfies at criticality the same recursion as in Eq. [34] 

^^f^+f^' (36) 



The recurrence of Eq [M] coincides with the transfer matrix 

K 



z L+ i = Y, G ~ tiZ L ( 37 ) 



i=l 



for the partition function Zl of a directed polymer on a Cayley tree of branching ratio K = 2 with random bond 
energies ej. [23|,[24j]. The differences with Eq. [34] we are interested in are the following : 

(i) the partition function Zl in Eq [37] is positive by definition, whereas here the random perturbation 8 n in Eq. 
1341 are a priori of arbitrary sign. Eq. [M] is thus more related to the case of a directed polymer model with complex 
weights studied in [25| . 

(ii) the weights e~ l3ei associated to the bond energies e, in Eq. [37jare now random weights % distributed with the 
fixed point distribution Q c (z) 

-> ^ (38) 

In particular, these weights present a broad power-law tail in \j z\ + ^ in contrast with the usual case where the energies 
ei are Gaussian. 

The difference (ii) turns out to be very important as we now explain. 
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2. Tails analysis 



Let us consider the first iteration 



Si 



(i) 

^(2)^^5(1) 



y (2) 
b "° 



(39) 



Suppose we start with a narrow distribution Vq(Sq) for the random initial perturbation So- The distribution "Pi((5i) 
after the first iteration will nevertheless present power-law tails inherited from the fixed point distribution Q c {z c ) ~ 
<I>(lnz c )/z^ +A1 with < /i < 1 (Eqs l25l and l33|) . More precisely, the tail in the region Si — > +oo is dominated by the 
events where zi 1 ^ is large with 5g 2 ^ > or where z£ ' is large with S^ > 0, and one obtains 



Vi(Si) ~ 2 / dz c Q c (z c ) / d5 Po(5o)5 

5i — >+oo 



+ OC 



<5i — -So 
b 



$(lnJi 



2 



dS Vo(So) <5£ 



Similarly, the left tail reads 



Vi(Si 



It is then clear that by iteration all distributions V n {S n ) will present these power-law tails 
Since we are looking for the Lyapunov exponent v governing the typical growth of the perturbation 

Sn+l v 



(40) 
(41) 

(42) 

(43) 
(44) 



it is convenient to rescale the iteration of Eq. [M] by the factor e v 



Vn+i = e 



,(2) 



..(2) 



(45) 



and to ask that the probability distribution P n (y) converges as n — > oo towards a stable distribution P x (y) presenting 
the tails (Eq[ 



y— *-±oo 



B 



(46) 



Reasoning as before, a large value of y n +i corresponds to a large value of one of the four variables (yn \ Dn \ z^), 
and one obtains the following equations 



B- 



B 



zr 



2e 






r + co 




bf 1 




Jo 


2e 


-LLV „M 


2e-^ 


f, 






h 6^ 


J — OO 



^Poo(y) y" 
^oo(y)|yr 



(47) 
(48) 



This shows that positive perturbations (initial distribution Vo(So < 0) = 0) or symmetric perturbations (symmetric 
initial distribution Vo(Sq) = Po{—S )) actually lead to the same Lyapunov exponent v 



2z£ 2 
bf* + 



~ — + — I dy ~B — V 



1+ ¥ L dy ^r y 



(49) 
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where we have used Eq. [28] The first term corresponds to the usual term for the velocity of the travelling wave 
approach [23], [24|, whereas the second term originates from the broad distribution of weights (z c /b). Its physical 
meaning is the following : in the usual case of a narrow distribution of the weights, the travelling wave approach 
allows to compute the velocity in terms of the weight statistics alone, because one can write a closed equation for 
the tail of the process [H, [24|; in the present case, the tail of the process does not satisfy a closed equation, because 
the broadness of the weight distribution induces some interaction between the tail and the bulk of the process : the 
second term in Eq. 29] represents the influence of the bulk of the distribution (y) onto the tail of exponent /j,. 

The exponent v describing the power-law growth 5 n ~ £n ~ e™ reads in terms of the Lyapunov exponent 

(50) 

v 

Note that the presence of the second term in Eq. 05] is crucial to obtain a finite exponent v : without this second 
term, the Lyapunov exponent v would vanish (v = 0) and the correlation length exponent would diverge {y — oo). 

III. NUMERICAL STUDY OF THE WETTING TRANSITION 
A. Numerical method 

We have performed numerical simulations with the so-called 'pool-method' which is very much used for disordered 
systems on hierarchical lattices [T3, EH : the idea is to represent the probability distribution P n {F n ) of the free-energy 
F n = —Tin Z n at generation n, by a pool of N values {F^, ..,Fn N ^}- The pool at generation (n+ 1) is then obtained 
as follows : each value F^^ is obtained by choosing two values at random from the pool of generation n and by 
applying the renormalization Eq. [7] 

The results presented in this Section have been obtained for the branching ratio b = 5, with a pool number 
N = 4.10 7 , with initial Gaussian energies (Eq. [3]). The corresponding annealed temperature is 

T ann = 1 ~ 0.60055 (51) 

Finally, the relation between the true free-energy F n = —ThxZ n and the reduced free-energy /„ = — T\nz n used 
in the previous section is simply (Eq[5]) 

F n = f n -T\nY n (52) 

where Y n = fe L " _1 does not contain any disorder. As a consequence, the two free-energy distributions have the same 
width AF n — Af n , and the same tail properties. 

B. Flow of the free-energy width AFl 

The flow of the free-energy width AFl as L grows is shown on Fig. [2] for many temperatures. One clearly sees the 
two attractive fixed points on this log-log plot. 

For T > T c , the free-energy width decays asymptotically with the exponent uj' w (b) introduced in Eq. [17] 

/ £ \-"w( b =5) \n(b 2 /2) 

AF(L) ~ with uj' w (b = 5) = -A^±J- = 1.8219.. (53) 

where £,p{T) is the corresponding correlation length that diverges as T — > T+. 

For T < T c , the free-energy width grows asymptotically with the exponent 1/2 (see Eq. [20l ) 

where £p (T) is the corresponding correlation length that diverges as T — > T~ . 

The critical temperature obtained by this pool method depends of the pool, i.e. of the discrete sampling with 
N values of the continuous probability distribution. It is expected to converge towards the thermodynamic critical 
temperature T c only in the limit N — * oo. Nevertheless, for each given pool, the flow of free-energy width allows a 
very precise determination of this pool-dependent critical temperature, for instance in the case considered 0.52415 < 
rppooi ^ 0.52416, which is significantly below the annealed temperature of Eq 1511 
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FIG. 2: (Color online) Wetting transition : log-log plot of the width AF(L) of the free-energy distribution as a function of L, 
for many temperatures. 



C. Divergence of the correlation lengths f F (T) 
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FIG. 3: (Color online) Wetting transition : Correlation length £ F (T) as measured from the behavior of the free-energy width 
(Eqs I53I and I54p (a) ln£^(T) as a function of T (b) lnf J(T) as a function of In \T C — T\ : the asymptotic slopes are of order 
j/~6.2 



The correlation lengths £,$(T) as measured from the free-energy width asymptotic behaviors above and below T c 
(Eqs IS"3l and l54l ) are shown on Fig. [3] (a). The plot in terms of the variable ln|T c — T\ shown on Fig. [3] (b) indicate 
a power-law divergence with the same exponent 



f^(r) m oc m \T-T C \-" with i/~6.2 



(55) 



D. Histogram of the free-energy 



The asymptotic probability distribution Hp of the rescaled free-energy 



F-F, 



Xf 



av(L) 



AF(L) 



(56) 



is shown on Fig [5] for three temperatures : 

(i) the distribution is Gaussian both for T > T c and T < T c as expected from Eqs [T5] and [TSl 
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FIG. 4: (Color online) Wetting transition: Log-plot of the asymptotic distribution II f of the rescaled free-energy Xf = af(l) 

in the low-temperature phase (here T = 0.25), in the high-temperature phase ( here T = 1) and at criticality (here T^°° l - 
0.524155) 



(ii) at criticality, one clearly see that a left-tail develops in the region / — > — oo with tail exponent rj c = 1 in 
agreement with Eg 1231 The corresponding power-law exponent of Eq[26]of the fixed-point distribution of Eq. [25] is of 
order 



H ~ 0.45 



(57) 



The measure is not very precise because one clearly sees on Fig. Q]that on top of this power-law, there exists oscillations 
reflecting the discrete nature of the renormalization. But this value is anyway in the expected interval of Eq [33l 



E. Flow of the energy width 
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FIG. 5: (Color online) Wetting transition : Flow of the widths AE(L) of the energy distribution as L grows (a) In AE(L) as a 
function of InL for many temperatures (b) Comparison of In AE(L) , In AS(L) and In AF(L) as a function of InL at criticality 
(T c pooi = 0.524155). 



The flow of the energy width AE(L) as L grows are shown on Fig. [5] for many temperatures. For T > T c , we find 
that the width decays asymptotically with the same exponent ui^(b) as the free-energy (Eq[ 



AE(L) 



with ^(6 = 5) 



ln(b 2 /2) 
2 In 2 



1.8219. 



(58) 
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For T < T c , this width grows asymptotically with the exponent 1/2 as the free-energy fEq [54]) 

AE(L) ~ (59) 

(60) 

Exactly at criticality, the free-energy AF(L) width converges towards a constant, whereas the energy width grows 
as a power-law (see Fig [5] b) 

AE(L) ~ V" with y c ~ 0.16 (61) 
This exponent is in agreement with the finite-size scaling relation y c = \ jv with v ~ 6.2 (see Eqs l55[) 



IV. DIRECTED POLYMER MODEL ON DIAMOND HIERARCHICAL LATTICES 

In this Section we study the directed polymer model, whose partition function satisfies the quadratic renormalization 
of Eq [5] In contrast with the wetting case described above, the transition only exists in the presence of disorder. 
Since we are interested into the asymptotic distribution of the free-energy, it is convenient to rewrite the free-energy 
of a sample (a) of generation n as 

=lnZW =F~+A n u a (62) 
where u a is a random variable of zero mean and width unity 

ul=l (63) 

So A n represents the width 

A^^-iKf) 1 ' 2 (64) 



A. Low-temperature RG flow 



In the low-temperature phase, the width A„ of the free-energy distribution grows with n. So at large scale, the 
recursion is dominated by the maximal term in Eq. [5] 

F n+1 ~ min x < a < b (f^-V + Fi 2a >) (65) 

This effective low-temperature recursion coincides with the recursion of the energy Eq of the ground state studied in 

El 

[HI]. The whole low-temperature phase is thus described by the zero-temperature fixed-point. In particular, the 
width of the free-energy distribution grows as 

A n ~L£°W (66) 
where luq (b) is the exponent governing the width of the ground-state energy AE - LT {b) studied in [ljj. 



B. High-temperature RG flow 

In the high-temperature phase, the width A n of the free-energy distribution is expected to decay to zero. The 
linearization in A„ of the recursion of Eq. [5] yields 



(3F n . 



hi 



y2 e -p(Fi 2a - 1) +Fi 2a) ) 



a=l 



= 2[3F n - In 



(U2a-1 + U2a) + 0((3 2 A^) 



a=l 



2f3 ln - ln(6) + ^ ^2(u 2a ^ + u 2a ) + 0((3 2 A 2 n ) 



(67) 
(68) 
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The consistence with the scaling form of Eq. [62] at generation (n + 1) 



F n+1 = F n+ i + A n+1 u (69) 

yields 



F n+1 = 2F n -T\n{b) 
b 



A„+im = 2j(u2o-i + u 2a ) (70) 



The normalization condition of Eq. [S3] yields 



A n+ i = y^A n (71) 



and 



1 b 

u = -j=y2{u 2a -i + u 2a ) (72) 

v a— 1 

For b > 2, this high-temperature phase exists, the probability distribution of the free-energy converges to a Gaussian. 
The width decays as the following power-law of the length L n = 2™ 

A„ocL- w ~W with Woo (6) = ^-L (73) 
In this regime, the rescaled variable u evolves according to Eq. 1721 and thus becomes Gaussian upon iteration. 

C. Analysis of the critical point 

At criticality, to avoid the high-temperature and low-temperature described above, the width A„ should converge 
as n — > oo towards a finite value A c . In particular, the fluctuating part of the free-energy 

= F^ -Fn = A c u^ (74) 

should remain a random variable of order O(l), of zero mean, distributed with some scale-invariant probability 
distribution P c (f) defined on ] - co, +oo[. Equivalently, the variable 

4«) = e-M^ = e -/»oA e uW (75) 
should have a scale-invariant probability distribution Q c {z) defined on [0, +oo[, with 

r + co 

\nz= I dzlnzQ c (z) = (76) 

The recursion for these variables z n reads 



a=l 



where 



B= lim (F n+1 -2F n ) (78) 

n — >oo 

should be finite ( otherwise the recursion of Eq. 1771 would not lead to a non-trivial stationary distribution Q c (z)), 
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1. Left-tail behavior of the free- energy distribution 

Let us introduce the left-tail exponent r] c 

lnP c (/) ~ - 7 (_/)*> + ... (79) 

f—>—oo 

where (...) denote the subleading terms. 

In the region where / — > — oo, one has effectively the low-temperature recursion of Eq. [65] 

/^m i n 1 < a < b (/( 2a - 1 )+/( 2Q )) (80) 

As in the wetting case, a saddle-point analysis shows that the only stable left tail exponent is 

Vc = 1 (81) 

So the distribution P c {f) decays exponentially 

P c (f) , * e<i\...) (82) 

f—>—oo 

This means that the corresponding distribution Q c (z) of z = e~P"f presents a power-law tail 

z^+oo z 

with some exponent 



(84) 

and where <&(lnz) represents the subleading terms. 
The first moment is fixed by the recursion of Eq 1771 

B 

Z=f (85) 
As a consequence, the exponent fi of the power-law of Eq [83] should satisfy 

p > 1 (86) 
and the parameter B representing the correction to extensivity (Eq I78[) is determined by 



B = bz = b J dzzQ c (z) (87) 
in terms of the fixed point distribution Q c {z). 

2. Analysis in terms of multiplicative stochastic processes 

Again, as explained in section III D 21 it is instructive to analyse the recursion of Eq. [77] from the point of view 
of multiplicative stochastic processes (see Appendix [A]) . The condition of Eq. IA2I translates here into the following 
condition using Eq. [7S] 



In | = -lnB= -\n{bz) < (88) 
The condition of Eq. IA4I translates into the following condition for the exponent fj, introduced in Eq. [531 

~Ty _ 2 JP _ 2 J+°°dzQ c (2 



2b 



j x °° dzQ c (z)z 



The argument is similar to Eqs IA5llA6[ the additional factor of 2b coming from the fact that z large corresponds to 
one of the z^' being large. This condition means in particular that the subleading term $(z) in Eq. 1251 should ensure 
the convergence at (+oo) as in the wetting case (see Eqsl2T)landl3"Tj). 
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Reminder on transitions of integer moments 



Let us now briefly recall the behaviors of the first moments discussed in [l4[ • From Eq [5j one obtains [l4 



-'ra+l 



b(Z^) 2 (90) 



2 



Z 2 n+i =b[Zl) +b(b-l)(Z n ) 4 (91) 
so that the ratio of the moments of the rescaled variable z defined in Eq [75] 

s # = i§ (92) 

follows the closed recursion [l4| 

r 2 (n+l) = ^^ (93) 

that actually coincides with Eq. [5] for the pure wetting model. The repulsive fixed point r* = b — 1 allows to define 
the temperature T 2 via r 2 (rt = 0) = b — 1 [1J|: for T > T 2 , the ratio r2(n) flows to 1, whereas for T > T 2 , the ratio 
r 2 (n) flows to (+oo). Similarly, the RG flow of ratios corresponding to higher moments have been studied in [LH ]. 
with the conclusion that for generic b (more precisely b > 2.303...) their transition temperatures are higher than T 2 . 
Since we already know /i > 1 (Eq 155]). we have to distinguish two cases 

(i) if the tail exponent /i satisfies 1 < fi < 2, then the second moment diverges at criticality z 2 = +oo and we have 
the strict inequality T c < T 2 . 

(ii) if the tail exponent \x satisfies fi > 2, then the second moment is finite at criticality. The only possible finite 
stable value is for the ratio r 2 is r 2 = b — 1. The critical temperature then coincides with the transition temperature 
of the second moment T c — T 2 . 

The scenario (i) is the most plausible, since the possibility (ii) would require some 'fine tuning' in some sense : as 
explained in the introduction, the temperature only appears in the initial condition (Eq[4]) of the renormalization ; 
any initial temperature T > T c flows towards the high temperature fixed-point, any initial temperature T < T c flows 
towards the low temperature fixed-point, so that T c is defined as the only initial temperature from which the critical 
distribution Q c (z) is accessible. The critical distribution Q c (z) has to satisfy the self-consistent equation of Eq 1551 
to be stable. If (ii) were true, the distribution Q c (z) should in addition satisfy a second completely independent 
condition r 2 = b — 1 , which seems unlikely. 

In conclusion, we expect that the exponent /i introduced in Eq. 1831 satisfies 

1< n < 2 (94) 
This is in agreement with the numerical simulations presented below in section fVl 



4. Right tail behavior of the free-energy distribution 



Let us introduce the right-tail exponent r\ c 

lnP c (/) ~ -i.f< + ... (95) 

In the region where / — * +00, one has effectively the high-temperature recursion 

/c-^/W (96) 

i=l 

where all free-energies are large. A saddle-point analysis with the right tail of Eq. [Ml shows that the only stable right 
exponent r/' c should satisfy b = 2 I?c_1 i.e. 
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D. Critical exponent v 



To compute the critical exponent v, we consider a small perturbation in the fluctuating part of the free-energy of 

- [3J^ = -[3 C (F^ - iQ = -0A c u^ + 4 o) (98) 
where 6^ represent the random perturbations of zero mean 

K = (99) 
Equivalently, these variables S n represent the perturbation at linear order of the variables of Eq [75] 

z^ = e ~^fL a) = Z M + g(p) (100) 
The linearization of the recursion of Eq. 1771 around the fixed point, yields 

2b (i) 

a n+1 ~££l_*« (101) 
i=l 

where are distributed with the critical distribution Q c (z). 

As in the wet ting case, the recurrence of Eq 11011 coincides with the recurrence describing a directed polymer on a 
Cayley tree [HI, l24j . with the following differences 

(i) the variables S n are random variables of zero mean fEq I99p . which is equivalent to a directed polymer model 
with random signs studied in [25| 

(ii) more importantly, the random weights are distributed with the fixed point distribution Q c (z) presenting a 
broad power-law tail in l/zl +fJl with 1 < fi < 2 (instead of < \l < 1 for the wetting case). 

Reasoning as in the wetting case, any narrow symmetric distribution Vq(6q) = Vq{—6q) will lead to power-law tails 
of index (1 +fi) after one iteration. The study of the evolution of these tails by iteration yields that the corresponding 
Lyapunov exponent v will be determined by an equation similar to Eq. 35] 

where we have used Eq. [89] in terms of the stationary distribution Poo{y) of the rescaled process associated to Eq llOll 

2b (i) 

yn+l = e-*Y,^ry$ ( 103 ) 

i=i 

The correlation length exponent then reads v = — . As in the wetting case, the presence of the second term in Eq 
11021 is crucial to obtain a positive v and a finite v. 



V. NUMERICAL STUDY OF THE DIRECTED POLYMER TRANSITION 

As for the wetting transition (see Section [III A[) . we have used the 'pool-method' with a pool number N = 4.10 7 to 
study the transition of the hierarchical lattice of branching ratio b = 5 with initial Gaussian energies (Eq. [3]). The 
exact bounds on the critical temperature are [14| 

T (b) = — r ~ 0.557... < TJb) < T 2 (b) = r ~ 0.849.. (104) 

[2 In b] 5 [ln(6-l)] 5 

In [3], the phase transition has been studied numerically via the specific heat and the overlap. In this paper, we 
characterize the transition via the statistics of free-energy and energy. As in the wetting case, this allows to locate 
very precisely the pool-dependent critical temperature and to measure the divergence of the correlation length £ (T) 
above and below T c . 
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A. Flow of the free-energy width 

The flow of the free-energy width AFl as L grows is shown on Fig. [6] for many temperatures. One clearly sees 
the two attractive fixed points. For T > T c , the free-energy width decays asymptotically with the exponent uj^b) 
introduced in Eq. [73] 



AF(L) 



L 



\n(b/2] 

with uj^Vb = 5) = \' ' = 0.6609. 
y ' 21n2 



(105) 



where Cp(T) is the corresponding correlation length that diverges as T — > T+. 

For T < T c , the free-energy width grows asymptotically with the exponent Wo(6) of the ground-state energy 
distribution 



AF(L) 



L 



with u> (b = 5) ~ 0.186. 



(106) 



where (T) is the corresponding correlation length that diverges asT-> T~ . 

For each given pool, the flow of free-energy width allows a very precise determination of the pool-dependent critical 
temperature, for instance in the case considered 0.77662 < Tf° o1 < 0.77666 which is significantly below the upper 
bound T 2 of EqUDS 



In AF 




FIG. 6: (Color online) Directed polymer transition : log-log plot of the width AF(L) of the free-energy distribution as a 
function of L, for many temperatures. 



B. 



Divergence of the correlation lengths £ F (T) 



The correlation lengths £p(T) as measured from the free-energy width asymptotic behaviors above and below T c 
(Eqs 11051 and 11061 ) are shown on Fig. [7](a). The plot in terms of the variable In \T C — T\ shown on Fig. 0(b) indicate 
a power-law divergence with the same exponent 



f!(T) m oc m \T-T C \-" with i/~3.4 



(107) 



Histogram of the free-energy 



The asymptotic probability distribution 11^ of the rescaled free-energy 



F-F 



Xp 



av(L) 



AF(L) 



(108) 
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FIG. 7: (Color online) Directed polymer transition : Correlation length $,p(T) as measured from the behavior of the free-energy 
width (Eqs ll05l and [T06ll (a) ln£^(T) as a function of T (b) ln£^(T) as a function of In \T C — T\ : the asymptotic slopes are of 
order v ~ 3.4 




FIG. 8: (Color online) Directed polymer transition : Asymptotic distribution n> of the rescaled free-energy xf = af(l)' > m 
the low-temperature phase (here T = 0.4), in the high-temperature phase ( here T — 2) and at criticality (here T^°° l — 0.77665) 
(a) Bulk representation (b) Log-representation to see the tails. 

is shown on Fig [8] for three temperatures : 

(i) for T > T c , it is a Gaussian in agreement with Eq 1721 

(ii) for T < T c , it coincides with the ground state energy distribution. 

(iii) at criticality, one clearly see that a left-tail develops with tail exponent rj c = 1 in agreement with Eg 1811 The 
corresponding power-law exponent of Eq [25] of the fixed-point distribution of Eq. [2H is of order 

(109) 

Again this measure is not precise as a consequence of the unknown logarithmic correction in Eq. 1251 but it is in the 
expected interval of EqIMl 

D. Flow of the energy and entropy widths 

The flow of the energy width AE(L) as L grows is shown on Fig. [5] (a) for many temperatures ( the flow of the 
entropy width AS(L) is very similar at large scale). For T > T c , we find that these widths decay asymptotically with 
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the same exponent u>oo(b) as the free-energy fEq I105[) 

AE{L) ~ L-^ {b) (110) 
AS(L) ~ L~ u °°^ (111) 

For T < T c , in agreement with the Fisher-Huse droplet scaling theory for directed polymers (26|, we find that these 
widths grow asymptotically with the exponent 1/2 which is bigger than the free-energy exponent u>o(b) (Eq. 1 1 06|) 



AE{L)~L* (112) 
AS(L)~L^ (113) 

Exactly at criticality, the free-energy AF(L) width converges towards a constant, whereas the energy and entropy 
widths grow as power-laws (see Fig [9] b) 

AE(L) ~ L Vc ~ AS(L) with y c - 0.29 (114) 

This exponent is in agreement with the finite-size scaling relation y c = l/i> with v ~ 3.4 (see Eqs ll07|) 

E. Divergence of the correlation lengths ^(T), £g(T) 

According to the Fisher-Huse droplet scaling theory of spin-glasses |26[, the singularities of the widths of energy 
and entropy as T ->■ T" is given by (L/^(T)) 1 / 2 /(T c - T). We thus define the correlation lengths £e( T ) and £j (T) 
by 

A£(£) ~ - J— 2 (115) 



AS(L) ~ — ^— ( — ^— ^) 2 (116) 
Similarly, for T > T c , we define the corresponding correlation lengths £e(T) and by the equations 



-aJoo(fe) 

2 / r \ -"oo (6) 



A£(L) ~ -J— f-I^— ] 'I ITi 



^'-r^lifm 1 (118) 
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FIG. 10: (Color online) Directed polymer transition : Correlation length Ce(T) (circles below and square above) and £f (T) 
(triangles below and diamond above) as measured from the behavior of the energy and entropy widths (a) ln^_e(T) and ln£s(T) 
as a function of T (b) ln^_B(T) and ln^s(T) as a function of In \T C — T\: the asymptotic slopes are of order v ~ 3.4 as in Fig. [7] 



The correlations lengths are shown on Fig. [10] (a) The plot in terms of the variable In \T C — T\ shown on Fig. \W\ 
(b) indicate a power-law divergence with the same exponent as in Eq. 11071 

^(T) r oc r \T-T C \-" with ^-3.4 (119) 



Histogram of the energy 



The asymptotic probability distribution He of the rescaled energy 



E-E, 



x E 



av(L) 



AE(L) 



(120) 



is shown for three temperatures on Fig 1111 

(i) outside criticality, both for T > T c and T < T c , these distributions are Gaussian. 

(ii) at criticality, the distribution is strongly non-gaussian and asymmetric, with a left-tail of tail exponent ?y c = 1 




in n 




-14 -12 -10 



E — E a v 



in the 



FIG. 11: (Color online) Directed polymer transition : Asymptotic distribution lie of the rescaled energy x 
low-temperature phase (here T = 0.4), in the high-temperature phase (here T — 2) and at criticality {T^°° l = 0.77665) (a) 
Bulk representation (b) Log-representation to see the tails 
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VI. COMPARISON WITH CORRESPONDING RESULTS ON HYPERCUBIC LATTICES 



Since the exact renormalizations on the diamond lattice can also be considered as approximate Migdal-Kadanoff 
renormalizations for hypercubic lattices, it is interesting to discuss whether the results obtained for the wetting and 
the directed polymer on the diamond lattice are qualitatively similar to the results for hypercubic lattices. 



A. Similarities for T < T c 



The whole low-temperature phase of disordered systems is usually characterized by the zero-temperature fixed point 
where disorder fluctuations dominate. For the disordered polymer models considered in this paper, the free-energy 
fluctuations grow as power-law of the length both for the diamond lattice and for hypercubic lattice 

AF(L, T < T c ) oc L^" (121) 

L — >oc 

where loq is the exponent governing the fluctuations of the ground state energy Eq(L). In the wetting case, this 
exponent has the simple value uj$ ett — 1/2 that reflects the normal fluctuations of the L random variables defining 
the random adsorbing energies along the wall. In the directed polymer case, the exponent loq is non-trivial because 
the ground state configuration is the result of an optimization problem. 



B. Differences for T > T c 



The high temperature phase of disordered systems is characterized by bounded disorder fluctuations, but these 
fluctuations are not of the same order on diamond lattices and on hypercubic lattices. More precisely, for the disordered 
polymer models considered in this paper, the free-energy fluctuations decays as a power-law on the diamond lattices, 
whereas they remain of order O(l) on hypercubic lattices 

Diamond: AF(L,T>T C ) cx L^^ h) (122) 
Hypercubic: AF(L,T > T c ) cx 0(1) (123) 

L — ►co 

This difference seems to come from the boundary conditions : (i) on the diamond lattice, the polymer is fixed at the 
two extreme points, but by the iterative construction of the lattice, the coordinence of these two extreme points grow 
with the number n of generations, so that it is possible to have a very efficient averaging even near the boundaries (ii) 
on the hypercubic lattices, the boundary conditions are sufficient to produce free-energies fluctuations of order 0(1) 
: the fixed origin has a finite coordinence, and the fluctuations of order O(l) of the random variables near this origin 
do not disappear as L — > oo. 



C. Differences at criticality 



On the hierarchical lattice, the free-energy fluctuations of the disordered polymer considered here are of order O(l) 
at criticality 

Diamond: AF(L,T = T C ) oc 0(1) (124) 

L — >oo 

and it is the only possibility in the presence of exact renormalizations : if the free-energy width is growing, the flow 
will be attracted at large scale towards the zero-temperature fixed point of Eq. 11211 whereas if free-energy width is 
decaying, the flow will be attracted towards the high-temperature fixed point of Eq. 11231 On hypercubic lattices, the 
free-energy fluctuations AF(L, T c ) ~ L Wc at criticality are expected to be governed by a vanishing exponent uj c = 0, 
but they are not necessarily of order 0(1) because logarithms cannot be excluded, and have actually been found 
for the directed polymer transition as we now explain. Forrest and Tang (27j have conjectured from their numerical 
results on a growth model in the KPZ universality class and from the exact solution of another growth model that 
the fluctuations of the height of the interface were logarithmic at criticality. For the directed polymer model, this 
translates into a logarithmic behavior of the free energy fluctuations at T c 

Hypercubic: AF^ d P (L, T = T c ) oc (lni) CT (125) 
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where the exponent was measured to be a = 1/2 in d = 3 [U [H, Further theoretical arguments in favour of 
this logarithmic behavior can be found in (30l . l3l| . So the scaling of free-energy fluctuations at criticality seem to be 
different on hypercubic lattices and on diamond lattices. 

Another related issue concerns the location of the critical temperature T c with respect to upper bound T 2 ■ 

(i) on the diamond lattice, the ratio r 2 of Eq[92]is finite at T 2 , the ratio z = Z / Z ann is a finite random variable at 
T c , but the probability distribution of the corresponding partition function presents a power law tail of index (1 + fi) 
with 1 < \l < 2 (Eq. [M]) , leading to the strict inequality T c < T 2 

(ii) on hypercubic lattices, the location of T c with respect to T 2 is still controversial. In [32j], we have argued that 
T c = T 2 in dimension d — 3, because the divergence of r 2 ~ e alnL at T 2 is compatible with the logarithmic free- 
energy fluctuations of Eq. 11251 provided the rescaled distribution of free-energy involves a left-tail exponent r\ c > 1 , 
as measured numerically in 29]. And in 33], we have found clear numerical evidence from the statistics of inverse 
participation ratios that the derealization transition takes place at T 2 . However, other arguments are in favor of the 
strict inequality T c < T 2 in finite dimensions : a new upper bound T* < T 2 was proposed in 1 + 3 (34| , and in [35[ 
the location of T c with respect to T 2 was shown to depend upon dimension and probability distribution of the bond 
energies. In particular for the gaussian distribution , th e result T c < T 2 is obtained only for d > 5 [35[, but not for the 
case d = 3 considered in numerical simulations [29l . |33| . 

For the wetting transition in 1 + 1 dimension, we are not aware of results concerning the scale of free-energy 
fluctuations at criticality. 

This comparison between the diamond lattice and hypercubic lattice can be summarized as follows. The free- 
energy fluctuations present analogous power-law behaviors in the low-temperature phase (Eq. I12ip but have different 
behaviors in the high temperature phase fEq |123p . At criticality, the free-energy fluctuations could also scale differently 
if logarithmic contributions are present on regular lattices. 

VII. CONCLUSION 

In this paper, we have studied the wetting transition and the directed polymer derealization transition on diamond 
hierarchical lattices. These two phase transitions with frozen disorder correspond to the critical points of quadratic 
renormalizations of the partition function. We have first explained why the comparison with multiplicative stochastic 
processes allows to understand the presence of a power-law tail in the fixed point distribution P c {z) ~ $(z)/z 1+M 
as z — > +00 ( up to some sub-leading logarithmic function $(z)) so that all moments z n with n > fj, diverge. The 
exponent fj, is in the range < \i < 1 for the wetting transition ( the first moment diverges z = +00 and the 
critical temperature is strictly below the annealed temperature T c < T ann ) and is the range 1 < ji < 2 for the 
directed polymer transition ( the second moment diverges z 2 = +00 and the critical temperature is strictly below the 
transition temperature T 2 of the second moment.) We have then obtained that the linearized renormalization around 
the critical point, which determines the exponent ^, coincides with the transfer matrix describing a directed polymer 
on the Cayley tree, where the random weights determined by the fixed point distribution P c (z) are broadly distributed. 
We have shown that it induces important differences with respect to the usual travelling wave solutions concerning 
more narrow distributions of the weights [23|, 0, HE}, where the selected velocity only depends on the tail region. 
Note that travelling waves also appear in other renormalization approaches of random systems [36| . Finally, we have 
presented detailed numerical results on the statistics of the free-energy and of the energy as a function of temperature 
for the wetting and the directed polymer transition on the diamond hierarchical lattice with branching ratio 6 = 5. 
In particular, we have shown that the measure of the free-energy width AF(L) yields a very clear signature of the 
transition and allows to measure the divergence of the correlation length ^(T) both below and above T c : (i) for 
T < T c , the free-energy width is governed by the zero-temperature exponent ujq via AF(L) ~ (L/£_(T)) W °; (ii) for 
T > T c , the free-energy width is governed by the high-temperature exponent via AF(L) ~ (L/^ + (T)) _ "°° . From 
the point of view of histograms, the development of a left tail with exponent rj c = 1 at criticality is very clear and 
different from histograms with exponent r\ > 1 outside criticality. 

APPENDIX A: REMINDER ON MULTIPLICATIVE STOCHASTIC PROCESSES 

Multiplicative stochastic processes appears in many contexts, in particular in one-dimensional disordered systems, 
such as random walk in random potentials [3?i . l38l . I39I ] or random spin chains [4p| . l4lj ] In this Appendix, we recall some 
useful results concerning the following recurrence of random variables X n 



X n +1 — UnXn + b n 



(Al) 
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where (<z n , b n ) are positive independent random numbers. The condition to have a stationary probability distribution 
P 00 (X)is 

hT^ < (A2) 
The most important property of P^X) is that it presents a power-law tail 

where the exponent /i > is determined by the condition [53, HH, H^, H3, 5l[ 

= 1 (A4) 
To understand where this condition comes from, one needs to write that P aD (X) is stable via the iteration of Eq. 

ED 

P o{X)= [ daV{a) [ db^(b) [ dYP oc (Y)5 (X - (aY + b)) = [ daV(a) [ dbip(b)^^^- (A5) 



where V(a) and i/)(b) are the probability distributions of a n and b n respectively. The stability of the power-law tail 
of EqEHlin the region X — > +oo yields at leading order 

( r daV(a) [dty(b)a*^-=W (A6) 



x l +^ J J X 1+ ^ X 1+ v 

yielding the condition of Eq. IA4I 
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